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Abstract 

Spectral and timing studies of Suzaku ToO observations of two SGRs, 1900+14 and 1806—20, are pre- 
sented. The X-ray quiescent emission spectra were well fitted by a two blackbody function or a blackbody 
plus a power law model. The non-thermal hard component discovered by INTEGRAL was detected by the 
PIN diodes and its spectrum was reproduced by the power law model reported by INTEGRAL. The XIS 
detected periodicity P = 5.1998±0.0002s for SGR1900+14 and P = 7.6022 ± 0.0007s for SGR1806-20. 
The pulsed fraction was related to the burst activity for SGR 1900+14. 
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1. Introduction 

Magnetars (neutron stars with a super strong surface 
magnetic field B ~ 10 15 G - e.g., Duncan & Thompson 
1992; Paczyhski 1992) have received considerable atten- 
tion recently. The magnetar's field strength exceeds the 
critical field B c = mlc 3 /ch « 4.4 x 10 13 G, and hence the 
non-linear effects of quantum electrodynamics must be 
considered in the processes of radiation transfer. Several 
studies have now found X-ray sources which are magnc- 
tar candidates, namely the soft gamma repeaters (SGRs) 
and anomalous X-ray pulsars (AXPs). They display X-ray 
quiescent emission with luminosities L = 10 33 -10 35 ergs -1 
(e.g., Woods & Thompson 2006). The SGRs and a 
few AXPs also exhibit "short bursts" with typical du- 
rations At ~ 100 ms and super-Eddington luminosities 
L ~ lO^ergs" 1 . One of the most remarkable burst phe- 
nomena is giant flares from SGRs. Three SGRs have emit- 
ted giant flares over the past three decades. The most en- 
ergetic one had L~ 5 x 10 47 crgs-\ from SGR1806-20 
on 27 December 2004 (Terasawa et al. 2005). The giant 
flares have been proposed as candidates for some short 
gamma-ray bursts (GRB) with durations At ^ 2s (e.g., 
Hurley et al. 2005). Delayed X-ray and radio emission 
following giant flares or short bursts display similar be- 
havior to GRB afterglows (Frail et al. 1999; Feroci et al. 
2003; Cameron et al. 2005; Nakagawa et al. 2008), which 
may also support this hypothesis. 



Two spectral models, a two blackbody function (2BB) 
and a blackbody plus a power law model (BB+PL), have 
been suggested for the quiescent X-ray emission of the 
magnctar candidates. A comprehensive study of many 
short SGR bursts detected by the High Energy Transient 
Explorer 2 (HETE-2) has shown that 2BB is preferable, 
although it is an empirical model (Nakagawa et al. 2007a). 
Although a physical model for the energy spectra remains 
controversial despite observations by several satellites, 
strong linear correlations between the lower and higher 
blackbody temperatures (/cTlt, fcTki), and between the 
lower and higher blackbody radii (i?L X , -Rht) were found, 
irrespective of the activity states (i.e., the X-ray quiescent 
emission or the short bursts), using 2BB (Nakagawa et al. 
2007b). There are various theoretical models which re- 
produce the spectra of the quiescent X-ray emission (e.g., 
Perna et al. 2001; Guver, Ozel & Lyutikov 2006). Even 
more exotic models have been proposed, such as that in- 
voking a p-star (Cea 2006). 

One piece of evidence in favor of the magnetar model 
is that magnetic pressure confines the photon-pair plasma 
from the giant flare in a (rather) small volume close to the 
neutron star (Thompson & Duncan 1995). One compli- 
cation is the fact that there has not been a secure direct 
measurement of the magnetic field of a magnetar candi- 
date. One possible such measurement is an absorption 
feature interpreted to be due to proton cyclotron reso- 
nance scattering which has been seen in a spectrum of the 
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precursor emission of a ~ 1.5s burst from SGR1806— 20 
(Ibrahim et al. 2002). This kind of feature can also be 
seen in the spectrum of X-ray quiescent emission from 
AXP1RXSJ170849-400910 (Rea et al. 2003). The lack 
of this feature in most magnet ar candidates may be ex- 
plained by a model which smears the absorption feature 
due to multiple cyclotron resonant scatterings in a stel- 
lar atmosphere and its magnetosphere (Giiver, Ozel & 
Lyutikov 2006). 

Several studies have reported periods P = 5-13 s and 
period derivatives P= 10~ n -10~ 13 s s" 1 for the magnetar 
candidates (Woods & Thompson 2006). It is known that 
pulse characteristics such as the pulsed fraction and pulse 
shape vary in time (e.g., Woods et al. 2007), and involve 
complex pulse morphology. However, little has been done 
to elucidate the underlying physics of the pulse character- 
istics. A recent theoretical study suggests that the pulse 
characteristics may be explained by trapped fireballs pro- 
duced by starquakes (Jia et al. 2008). 

This paper reports the broadband spectroscopy of two 
SGRs, 1900+14 and 1806-20, observed by Suzaku. The 
relation between pulsed fraction and burst rate will also 
be discussed. The distances are assumed to be 14.5 kpc for 
SGR 1900+14 (Hurley et al. 1999; Vrba et al. 2000) and 
8.7 kpc for SGR 1806-20 (Corbel et al. 1997; Corbel & 
Eikenberry 2004; Cameron ct al. 2005; McClure-Grifhths 
& Gaensler 2005; Bibby et al. 2008). 

2. Observations 

Suzaku (Mitsuda et al. 2007) has the capability of giv- 
ing high quality broadband spectra of astrophysical X- 
ray and gamma-ray sources, because of the great sen- 
sitivity of the X-ray imaging spectrometer (XIS; 0.2- 
12keV; Koyama et al. 2007) and an extremely low back- 
ground level of the hard X-ray detector (HXD; 10-700 keV; 
Takahashi et al. 2007). The XIS consists of one back- 
illuminated (BI) CCD and three front-illuminated (FI) 
CCDs. The HXD consists of PIN diodes (10-70 keV) 
and GSO units (40-600 keV). The XIS has an effective 
area of 370 cm 2 (BI) and 330 cm 2 (FI) at 8kcV. Although 
the effective areas are comparable to the EPIC on-board 
XMM-Newton, XIS has better energy resolution below 
IkeV: -50 cV (BI) and ~40eV (FI) at 0.525 keV. The 
HXD PIN also has better energy resolution, ^3keV, 
compared with FREGATE 1 on-board HETE-2 (5keV 
at 20kcV) and IBIS on-board INTEGRAL (7keV at 
100 keV; Winkler 1999). Suzaku can observe objects in 
the energy range 0.2-600 keV which covers a lower X- 
ray band compared with HETE-2 (2-400 keV; Ricker et 
al. 2003) and INTEGRAL IBIS (15kcV-10MeV; Winkler 
1999). Since the HXD field of view is narrow (34' x 34'), 
it has good sensitivity to weak sources. The sensitiv- 
ity of the HXD is 3 x 10~ 6 photons s" 1 keV" 1 cm" 2 at 
20keV (Takahashi et al. 2007), which is better than 
HETE-2 FREGATE (~ 10" 3 photons s" 1 keV" 1 cm" 2 at 

1 See http://space.mit.edu/HETE/fregate.html. 



100 keV; Atteia et al. 2003) and INTEGRAL IBIS 2 
(10~ 5 photons s^kc-V" 1 cm" 2 at 20kcV). 

When the SGRs 1900+14 and 1806-20 displayed in- 
tense bursting activity, scheduled target of opportunity 
(ToO) observations by Suzaku were performed. Table 1 
shows a summary of the observations. 

SGR 1900+14 began a very high burst activity on 29 
March 2006, and Swift recorded more than 40 bursts on 
that day (Israel et al. 2008). The Suzaku ToO observation 
started at 8:43 and ended at 21:52 on 1 April 2006 (UT). 
The X-ray quiescent emission was detected by the XIS, 
but there was no significant emission in the HXD data. 
No burst was found in the XIS and HXD light curves. 
There were no Swift bursts during the Suzaku observation 
(the bottom right panel in figure 3). 

A bright long burst from SGR 1806— 20 was detected 
by Konus-Wind and Suzaku-WAM on 26 March 2007 
(Golcnctskii et al. 2007). The Suzaku ToO observation of 
this source started at 15:08 on 30 March 2007 and ended at 
1:30 31 March 2007 (UT). Both the X-ray quiescent emis- 
sion and the non-thermal hard X-ray emission discovered 
by INTEGRAL (Molkov et al. 2005) were detected by the 
XIS and the HXD. Two dim short bursts were detected by 
Suzaku. One of them was a short burst localized by Swift 
at 16:14:38, and the other was detected only by Suzaku 
at 17:34:00. Since these bursts are too weak to perform 
spectral analyses, we do not consider them further here. 

3. Analysis 

3.1. Data Reduction 

The reduction of both the XIS and HXD event data 
(v2.0) were made using HEAsoft 6.4.1 software. The lat- 
est calibration database (CALDB: 20080616) was applied 
to the unfiltered XIS event data using xispi (v2008-04-10). 
Then the new data were filtered using the basic criteria 3 
and with a grade selection "GRADE = (0,2-4,6)" using 
xselect (v2.4a). After that, hot and flickering pixels were 
removed using cleansis (vl.7). Telemetry saturated time 
intervals were estimated by xisgtigen (v2007-05-14), and 
the time intervals were removed using xselect. The sizes 
of the source and background regions are summarized in 
table 1. The background regions of the FI CCDs for the 
SGR 1806— 20 observation are selected near the edge op- 
posite to the source region (because the object fell on the 
edge), while those of the other CCDs were selected near 
both edges. Light curves and spectra were extracted from 
the clean XIS event data using xselect. The response ma- 
trix files were created by xisrmfgen (v2007-05-14) and the 
ancillary response function files were generated by xissi- 
marfgen (v2008-04-05). The spectra were binned to at 
least 50 counts in each spectral bin by grppha. 

The barycentric correction was applied to the cleaned 
XIS event data by use of aebarycen. Then 0.2-12 keV 
light curves with 1 s binning (the minimum time resolu- 
tion of the XIS 1/8 window mode) were extracted using 

2 See http://integral.esa.int/integ_payload_imager.html. 

3 The XIS basic criteria were derived from 
http: / / suzaku.gsfc.nasa.gov/docs / suzaku / analysis /abc. 
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the source and background regions for the spectra. To find 
the most reliable period, timing analyses were performed 
by powspec and efsearch. A folded light curve was made 
with the best period by efold. 

The latest calibration database (CALDB: 20080602) 
was applied to the unfiltered HXD event data using hxdpi 
and hxdgrade. The PIN and GSO event data were ex- 
tracted from the newly calibrated data with basic cri- 
teria 4 using xselect. The tuned background event data 
(bgd_d) were utilized to estimate the non X-ray back- 
ground (NXB). The cosmic X-ray background (CXB) 
was estimated based on the HEAO-1 result (Boldt 1987). 
Then the good time interval (GTI) files were generated by 
merging the PIN/GSO GTI extension and the bgd_d GTI 
extension using mgtime. Light curves and spectra were 
extracted using xselect. For the source spectra, the dead 
time corrections were applied using hxddtcor. The mod- 
eled background event data (bgd_d) were utilized in our 
analyses. The spectra were binned to at least 30 counts 
in each spectral bin by grppha. 

Since both SGRs 1900+14 and 1806-20 are located 
in the galactic plane, the galactic ridge emission (GRE) 
should be considered as well as NXB and CXB. It is not 
easy to estimate the contribution of the GRE, so we used 
the HESS J1804-216.BGD data (Bamba et al. 2007) to 
estimate both the CXB and GRE for the SGR 1806-20 
observation. Note that all PIN diodes were operated 
at 500 V during the HESS J1804— 216_BGD observation, 
while half of the PIN diodes (UNITID<7) were operated 
at 400 V during the SGR 1806-20 observation. Since si- 
multaneous analysis of these two observations with differ- 
ent modes is not available, we concentrate on the 500 V 
data (i.e., UNITID>7). 

3.2. Spectral Analysis 

Since a two blackbody function (2BB) and a blackbody 
plus a power law model (BB+PL) have been suggested 
for the quiescent SGR X-ray spectra, we used them for 
our spectral analyses. An absorption model was applied 
to both models. 

The XIS spectra of SGR 1900+14 are acceptably fitted 
by 2BB and BB+PL, for which the spectral parameters 
are summarized in table 3, and the spectra are shown in 
figure 1 (a). 

The non-thermal hard X-ray emission component dis- 
covered by INTEGRAL (Molkov et al. 2005) was detected 
by the PIN diodes for SGR 1806-20. The emission is sig- 
nificant even if we increase the background flux by a factor 
of 1.02 (corresponding to the current background fluctu- 
ation of the PIN diodes). Since there are not enough 
statistics to investigate the spectral shape of the non- 
thermal hard emission, we performed a spectral analy- 
sis using a power law model (HPL) with a fixed index of 
T = 1.6 as measured by INTEGRAL. The 10-70 keV flux 
was F= 2.9~!ii 2 x 10 -11 erg cm -2 s _1 , where the quoted er- 
rors are 68% confidence levels. The unabsorbed 1-200 keV 

4 The HXD basic criteria were taken from 
http: / /www. astro. isas.ac.jp/suzaku/analysis/7step_HXD_20080f 



flux was estimated to be F ~ 4.3 x 10 11 erg cm 2 s 1 
using the best-fit spectral parameters which is lower by 
a factor of ~ 3 than the flux observed by INTEGRAL 
(F ~ 1.3 x 10~ 10 crgcm- 2 s- 1 ; Molkov et al. 2005). Then 
we performed the broadband spectral analyses of the 
XIS and the PIN diodes for SGR 1806-20. The spec- 
tra are well described by 2BB plus HPL (2BB+HPL) or 
BB+PL plus HPL (BB+PL+HPL). The index and flux 
of the non-thermal hard emission are fixed to L = 1.6 and 
F = 2.9 x 10~ n crgcm~ 2 s" 1 , respectively The broadband 
spectra of SGR 1806— 20 are shown in figure 1 (b). 

3.3. Timing Analysis 

The pulse period of SGR 1900+14 was searched for 
using the XIS 0.2-12 keV light curves. Two candidates 
with nearly the same confidence level were found from an 
epoch-folding analysis: P = 5.1998 + 0.0002 s {\ 2 = 24.8) 
and P = 5.2043 ± 0.0003s ( X 2 = 24.7). The folded light 
curve with the first period shows a pulsation (the left panel 
in figure 2), while the other case does not display a clear 
pulsation (not shown in figure 2). In addition, the first 
pulse period is consistent with that determined by XMM- 
Newton on April 1 2006 (Mereghctti et al. 2006), the same 
date as that of the Suzaku observation. The most plausible 
pulse period is therefore P = 5.1998 ± 0.0002s. The pulse 
period derivative between the XMM-Newton observation 
in September 2005 (P = 5.198346±0.000003s; Mereghetti 
et al. 2006) and the Suzaku observation in April 2006 is 
estimated to be P= (8.7+1.2) x 10 -11 ss -1 . This is lower 
by a factor of 0.8 than the first measurement ~ 100 days 
prior to the giant flare on 27 August 1998 (Kouveliotou 
et al. 1999). Using the folded light curve with the best 
period, the pulsed fraction is found to be Pf = 16 + 3%. 
Here P f = (C max - C min )/ (Cmax + C min ), where C max and 
C m i n are the measured count rates at the maximum and 
at the minimum, respectively. 

The pulse period of SGR 1806— 20 from an epoch- 
folding analysis is found to be P = 7.6022 ± 0.0007s 
(X 2 = 53.1) using XIS 0.2-12 keV light curves. The folded 
light curve displays the pulsations, as shown in the left 
panel in figure 2. Comparing this with the reported 
pulse period derived from the ~XMM-Newton observation 
on September 10 2006 (P = 7.5891 ± 0.0002 s; Esposito 
et al. 2007), the pulse period derivative is estimated to 
be P = (7.5 + 0.4) x 10~ 10 ss _1 . This is consistent with 
previously reported values derived from observations after 
the giant flare on 27 December 2004 (e.g., Woods et al. 
2007). Using the folded light curve, the pulsed fraction is 
found to be P f = 8 ± 2%. 

4. Discussion and Conclusion 

The quiescent emission spectra of the SGRs 1900+14 
and 1806—20 measured by the XIS on-board Suzaku are 
well described by either a two blackbody function (2BB) 
or a blackbody plus a power law model (BB+PL). The 
2BB temperatures (/cTlt and A:Tht) and radii (i?£ T and 
_Bs^ T ) show good agreement with the /cTlt-^Tht and i? 2 T - 
-fi^HT correlations reported by Nakagawa et al. (2007b). 
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A shallow dip is apparent at around 2.3 keV in the 
SGR 1900+14 spectra as shown in figure 1 (a). However, 
the XIS hardware team has cautioned that the detector 
gain has some unresolved uncertainties around 2keV for 
the window mode 5 , so this dip might not be real. If we 
treat it as a line feature, however, it could be well rep- 
resented by a Gaussian shape absorption with E = 2.3 ± 
0.1 keV and a = 0.2_^ i keV, where E is the line energy and 
a is the line width. If the dip is due to the fundamental ab- 
sorption feature of proton cyclotron resonant scatterings, 
the dipole magnetic field is B = 3.7 x 10 14 (£'/2.3kcV) G. 
We conducted an absorption line search for both SGRs 
1900+14 and 1806—20 using a narrow Gaussian absorp- 
tion line model with the line center as a free parameter 
ranging from 0.2 keV to 12keV. No feature was found ex- 
cept for the above mentioned E ~ 2.3 keV shallow dip of 
SGR 1900+14. The upper limits of the feature are r < 0.3 
(SGR 1900+14) and r < 0.4 (SGR 1806-20). 

The non-thermal hard X-ray emission of SGR 1806—20, 
discovered by INTEGRAL in a 1.6 Ms long observation, 
was detected by the PIN diodes on-board Suzaku with a 
short exposure time (~ 20 ks) because of the good sensitiv- 
ity and the very low background level of this instrument. 

The periodicities of SGRs 1900+14 and 1806-20 (table 
2) are detected despite the short exposure time (~ 20 ks), 
again because of the low background level of the Suzaku 
XIS. It has been recognized that the pulsed fractions of 
SGRs are not stable, and that their pulse shapes are com- 
plicated. One example is the fact that the pulsed fraction 
displayed a clear drop-off around the day of the giant flare 
from SGR 1806-20 on 27 December 2004 (Woods et al. 
2007) . To perform a deeper study of the pulse character- 
istics of SGR 1900+14, we compared the pulsed fractions 
to the X-ray fluxes and burst rate. The pulsed fractions 
and the X-ray fluxes have been derived both from our 
work (tables 2 and 3) and from the literature (Mereghetti 
et al. 2006; Esposito et al. 2007; Israel et al. 2008), while 
the burst rate was estimated from the list on the IPN web 
site 6 . Panel (a) in figure 3 shows their time histories over 
10 years from 1997 to 2007, where A, B and C indicate the 
days of the giant flare on 27 August 1998, the intermediate 
flare on 18 April 2001 and the extraordinary burst activity 
with more than 40 bursts on 29 March 2006, respectively. 
The results of the Suzaku observation are plotted with star 
symbols. 

There seems to be drop-off of the pulsed fraction around 
the days of B and C. Panels (b) and (c) in figure 3 give 
an expanded view of panel (a) for days B and C, respec- 
tively. In panels (b) and (c), the pulsed fractions display 
a clear drop-off after both active days (i.e., B and C). 
No burst activity was found after the day B for 9 days, 
and then the quiet burst activity appeared. Since the day 
C, the burst activity had been gradually disappearing. 
No bursts were reported by the IPN during the Suzaku 
observation (the star symbol in figure 3). The known 



bursts before and after the Suzaku observation were de- 
tected at 12:19:51 on 29 March 2006 and at 00:23:48 on 5 
April 2006, respectively. On the contrary, SGR 1806-20 
displayed one burst during the Suzaku observation, and 
moderate, steady burst activity before and after it. This 
drop-off is found in response not only to the giant flare 
but also to the extraordinary burst activity. (This implies 
that the pulsed fraction is related to burst activity.) One 
possible explanation is to assume that magnetic confine- 
ment stores a vast amount of energy in the neutron star 
atmosphere or magnctosphcre through energy injections 
such as starquakes, and that the pulsed fraction is then 
increased. After that the stored energy may be released 
as bursts either through reconnection, or some sort of in- 
stability such as the thermonuclear explosions of X-ray 
bursts. A recent theoretical study reports that a variety 
of pulse morphologies can be caused by trapped fireballs 
released by starquakes (Jia et al. 2008). 

We conclude by pointing out that monitoring the time 
variations of the pulsed fraction and the burst rate is 
an important element in revealing the burst mechanism. 
This will be achieved by wide field detectors such as the 
Monitor of All-sky X-ray Image (MAXI) on-board the in- 
ternational space station (Matsuoka et al. 1997). 

We would like to thank K. Hurley for helpful comments 
and suggestions to improve our paper. We are also grate- 
ful to the XIS team for useful discussions on the timing 
analyses and the calibrations of the window mode. This 
work is supported in part by a special postdoctoral re- 
searchers program in RIKEN. Y.E.N, is supported by the 
JSPS Research Fellowships for Young Scientists. 



The reports are available on 

http: / /www. astro. isas.jaxa.jp/suzaku/proccss/caveats/cavoats_xrtxis06. html. 

The burst list is available on 

http: / /www. ssl.berkeley.edu/ipn3/sgrlist.txt. 
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Fig. 1. Quiescent emission spectra 
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of SGR 1900+14 (a) and SGR 1806—20 (b). 
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Fig. 2. Folded 0.2-12 keV light curves for SGR 1900+14 (a) and SGR 1806-20 (b). 



Table 1. A summary of Suzaku ToO observations. 



SGR SeqNum* Observation Date (MJD) W 

Start End (ks) 



los^ 



Regions^ 
FI 



BI 



1900+14 401022010 53826.363 53826.911 17.0 1,2 2?3 X 5.'0/2.'3 X 6.'7 
1806-20 401021010 54189.631 54190.063 19.6 0,1,3 2.'2x5.'3/2.'2x3.'7 

* Suzaku sequence number. 

t Net exposure time for each observation. 

J XIS sensor number. 

§ Extracted source/background regions for front-illuminated (FI) and back-illuminated (BI) CCDs. 



2f3 x 7f0/2.'3 x 3:2 
2:2 x 7:0/2:2 x 3:5 
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Fig. 3. Panel (a): time variations of the 2-10 keV flux, the pulsed fraction and the burst rate from 1997 to 2007, where A, B 
and C indicate the days of the giant flare on 27 August 1998, the intermediate flare on 18 April 2001 and the extraordinary burst 
activity with more than 40 bursts on 29 March 2006, respectively. Panels (b) and (c): expanded plots around the days of B and C, 
respectively. The Suzaku observation is plotted with star symbols. 



Table 2. A summary of timing analyses. 



SGR Epoch (MJD TDB) P* P f 

1900+14 53826.001 5.1998 ± 0.0002 8.7 ±1.2 

1806-20 54189.0009 7.6022 + 0.0007 75 ±4 

* Fulse periods in units ot s. 

' Pulse period derivatives in units of 10 — 11 ss _1 . 
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Table 3. A summary of spectral parameters. The XIS spectra of SGR 1900+14 are fitted with 2BB or BB+PL, while the XIS 
and PIN diodes spectra of SGR 1806-20 are fitted with 2BB+HPL or BB+PL+HPL. The spectral parameters of the non-ther- 
mal hard emission of SGR 1806— 20 are not shown in the table because the index and the flux were fixed to T = 1.6 and 
F = 2.9 X 10 — 11 ergem" 2 s" 1 (see subsection 3.2), respectively. 



SGR 


Model 


(10 22 cm" 


- 2 ) (keV) 


Rlt (Rbb)^ 
(km) 


(keV) 


-Rht* 
(km) 




F n 


X 2 (d.o.f.) 


1900+14 


2BB 


2.2 + 0.3 


0.49 + 0.06 




i 5 +u.3 


0.36'^ 




4.1 + 0.5 


115 (105) 


1806-20 


BB+PL 

2BB+HPL 

BB+PL+HPL 


q 9+0.9 
°- z -0.7 

'• S -0.9 

s.o±?:S 


29+ - 17 

u - oy: -0.09 

51 + 09 
U ' O1 -0.13 


1 9 q + 71.8 

"■«*_ 9.3 

3.li?5 


2.7±f;i 




2-s±g:S 


4.1 + 0.5 
9.9 + 1.6 
9.9 + 1.6 


106 (105) 
196 (192) 
196 (192) 



* Nn denotes the photoelectric absorption with 90% confidence level errors. 

t fcTLT, &7ht and A;Tbb denote blackbody temperatures with 90% confidence level errors. 

Rlt, Rut and -Rbb denote the emission radius with 90% confidence level errors. 
§ r denotes the power law index with 90% confidence level errors. 

II F denotes the 2-10 keV flux in units of 10 -12 erg cm -2 s -1 with 68% confidence level errors. 
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